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Abstract 

We study the vortex nucleation in a finite Bose-Einstein condensate. Using a set of non-local 
and chiral boundary conditions to solve the Schrodinger equation of non-interacting bosons in a 
rotating trap, we obtain a quantitative expression for the characteristic angular velocity for vortex 
nucleation in a condensate which is found to be 35% of the transverse harmonic trapping frequency. 
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I. INTRODUCTION 

A prominent feature of a superfluid is the way it behaves under rotation . In contrast to 
a normal fluid, that rotates like a rigid body at thermal equilibrium, the thermodynamically 
stable state of a superfluid at low enough frequency does not rotate. At higher frequencies, 
a finite amount of angular momentum appears in the form of vortex filaments at which the 
superfluid density vanishes. The circulation of the velocity field flow evaluated on a closed 
contour that encircles the vortex is quantized . This is a consequence of the existence 
of a macroscopic wavefunction, whose phase changes by an integer multiple of 27r around 



the vortex fllaments. Atomic Bose-Einstein condensates (BEC) 



provide a reference 



system where the superfluid behaviour can be studied in the weak-coupling regime. 



Idled m tne weak-coupimg : 



the existence 



A set of recent experiments has demonstrated 
of vortices in atomic condensates. The rotational frequency at which the flrst vortex is 
nucleated has been shown to depend on parameters such as the trap geometry (aspect 
ratio), the nature and the characteristic time of the stirring beam, the number of atoms 
in the trap and may therefore vary from one experiment to another. The ratio of this 
characteristic rotational frequency to the transverse harmonic trapping frequency varies in 
these experiments from 0.1 12| to 0.7 ^. 



Most of theoretical studies of vortex nucleation in atomic BEC which preceded these 



experiments 
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determine the characteristic nucleation frequency of the 
flrst vortex from the criterion that the vortex-state becomes the minimum of the thermo- 
dynamic free energy of the system evaluated in the co-rotating frame. In 21| the global as 



well as local stability of the vortex state is discussed within the framework of the Bogoli- 
ubov theory jsj. Using the Thomas- Fermi approximation, the thermodynamic characteristic 



frequency can be expressed 3, llH in terms of the other system parameters. These works 



have been reviewed in 1221. The characteristic rotational frequency of the vortex nucleation 
thus obtained is generally lower than the values observed experimentally 

The presence of a vortex in a static condensate is associated with an extra energy relative 
to the vortex free state. Using the Thomas-Fermi approximation it can be shown that this 
energy is the highest when the vortex is located at the center of the trap and decreases mono- 
tonic ally as a function of the distance between the center of the vortex and the center of the 
trap QHQ- 



Within Thomas-Fermi approximation (without any boundary correction) 
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it can be shown that this energy vanishes at the boundary of the system. Thus we can 
define the energy of a vortex-state as a function of the distance d, namely the separation 
between the vortex and the center of the trap and E^{d) has a maximum at d = 0. For an 
increasing rotational frequency this maximum is shifted from the center to the boundary of 
the system 0, and at rotational frequencies higher than the thermodynamic charac- 
teristic frequency this leads to a surface energy barrier to the nucleation of a vortex. This 
explains why a vortex cannot be nucleated in a trapped condensate even though the trap 
is rotated at the thermodynamic characteristic frequency. The characteristic frequency for 
the vortex nucleation can be determined if one knows under what condition and at what 
rotational frequency a state carrying finite angular momentum will be transferred from the 
surface to the bulk of the condensate by overcoming the surface energy barrier. Theoretical 
studies in this direction have been done in 

0, H Q, 
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They are 



based on an analysis of the co 
namely the surface modes 



lective excitations localized at the surface of the condensate, 
|. These modes appear as shape deformations that carry a 
finite angular momentum about the axis of rotation [3l and have no radial node. In a 
rotated condensate these surface modes are excited and this leads to the vortex nucleation. 
A generalization of the Landau criterion J, 
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allows to determine at what rotational 



frequency, the surface mode corresponding to a given angular momentum quantum number 
is excited and leads to the nucleation of a vortex. The characteristic nucleation determined 
in this way agrees with the value experimentally observed [2^. Vortices are also nucleated 



by exciting a particular surface mode through a controlled trap deformation jSl, llJ| and this 
problem has been theoretically studied in [^jQl- 

Thus far theoretical studies of the problem of vortex nucleation in a trapped condensate 
has been conducted in the framework of interacting bosons using various approximation 
schemes. It is also known that the stable state of a set of non-interacting bosons in an 
axisymmetric harmonic trap does not have finite angular momentum along the z-direction 
for a rotational frequency less than the trap frequency. A finite amount of interaction makes 
vortex states energetically feasible at a lower rotational frequency. In this article we study 
the nucleation of a vortex from the surface to the bulk of a system by solving the one- 
particle Schrodinger equation with a set of non-local and chiral boundary conditions. To 
that purpose we start by discussing in section II the role of boundary conditions in a many- 
body problem. In section III, we review the problem of a trapped boson rotating at a given 
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frequency in an infinite plane and briefly discuss the corresponding energy spectrum in a 
finite disc while applying Dirichlet boundary conditions. In section (IV) we introduce the 
chiral boundary conditions for rotating bosons in a disc and we show how the Hilbert space 
splits into bulk and edge states. We subsequently analyze the nucleation mechanism, i.e. 
the transfer of angular momentum form edge to bulk. The characteristic angular rotation 
is given in terms of the trap frequency at which the first and then successive vortices are 
nucleated in the bulk. The variation of the size of the bulk region with increasing rotation 
is discussed and its physical implication is pointed out. 



II. ROLE OF THE BOUNDARY CONDITIONS IN A MANY BODY PROBLEM 

We start with the following hamiltonian corresponding to interacting bosons of mass 

m 

N ^2 

1=1 i,j 

where Eo is the ground state energy. Except for some specific cases, one does not know how to 
diagonalize this hamiltonian. Therefore some approximation schemes must be defined whose 
purpose is to obtain an effective quadratic hamiltonian. We may consider, for instance, the 
Feynman description 



3q | that accounts for the excited states under the form 

(j){ri,--- ,rN) = F<po{ri,--- ^vn) (2) 

where 0o is the exact but unknown ground state wavefunction such that Hmb4>o = 0, and we 
assume that F = fi^fi)- The writing of F as a sum over one body terms is exact for the 
non-interacting problem. For the interacting case, such a decomposition assumes that the 
interaction is adiabatically switched on. The wavefunction ^(r) is obtained by minimizing 
the energy 



where d^r = dri ■ ■ ■ dr^^. It can be shown that 3^ Q the effective energy E can be written 

as 

E=-Po^ I drnr)V'f{r) (4) 

where po is the ground state density. To obtain the corresponding spectrum, we have to 
impose boundary conditions on /(r). Assuming translational invariance, Feynman obtained 



that E = ^[^^], where S{k) is the structure factor. In a trapped condensate the trans- 
lational invariance is broken by the presence of a confinement potential. The function /(r) 
is related to the order parameter \E'(r). The choices for boundary conditions on the order 
parameter \E^(r) is broad as it depends on the nature of the confining potential and on the 
effective one body term generated by the interaction. It might be thus possible to take into 
account at least partly the effect of interactions by solving a linear Schrodinger like equation 
under suitable choices of boundary conditions. However in the absence of a specific relation 
between such boundary conditions and the effective interactions, these choices are generally 
guided by the nature of the problem. 

In this paper we study the vortex nucleation in a confined geometry by solving a one- 
particle Schrodinger equation for the condensate wavefunction \E^(r) with a set of non-local 
and chiral boundary conditions. Such boundary conditions have been proposed in order to 



following considerations. We know 
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deal with such non-linear problems 37 1. These boundary conditions are motivated by the 
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that the dispersion 



relation for the surface excitations determines the characteristic rotational frequency of 
vortex nucleation. The proposed boundary conditions are designed in order to provide a 
clear distinction between the bulk and edge states of a two dimensional rotating boson gas. 
The angular momentum quantum numbers of edge states as we shall see are higher than 
those of bulk states. Vortices are then nucleated by transferring a state from the edge to 
the bulk Hilbert spaces. 

III. A ROTATING TWO DIMENSIONAL BOSON GAS 
A. Hamiltonian and energy spectrum for the infinite plane 

We consider the hamiltonian of a trapped boson in a two dimensional domain rotating 
with a uniform angular frequency Q : 

r>2 1 

H = — + -mwV^ - QL, (5) 
2m 2 ^ ^ 

We define the vector potential 

A/ = / X r (6) 
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where / = (0, 0, /). Corresponding pseudo-magnetic fields may be defined by 

= V X = 2/ (7) 
where / is either a; or so that the hamiltonian (jS} rewrites under the two equivalent forms 



or 



H = -^(p - mA^f + {uj- n)L, (8) 
2m 

H = ^(p - mAnf + ^m{uj^ - (9) 

For uj = Q, this hamiltonian is the Landau hamiltonian of a charged particle in a transverse 
magnetic field written in the symmetric gauge. For this special value the centrifugal force 
just offsets the confinement and hence the bosonic system becomes unstable. 

The eigenfunctions and the eigenvalues En,i of the hamiltonian for the infinite plane 
are characterized by two integer quantum numbers n, I where n G N, / G Z. We define 
buj = and 6q = The solutions of the corresponding Schrodinger equation are 

^n,; W = C„,;rl'le^'^e-'^ iF,{a, |/| + 1; b^r^) (10) 

where a = ^^"'^"^ — | = —n with e = iFi{a,c;x) is the confluent hypergeometric 

function [s^l^] and Cn,i is a normalization constant. The eigenenergies are 

E^^i = hw{2n+\l\--l+l) (11) 
u 

For / < 0, the eigenvalues increase with increasing Q, whereas for / > 0, they decrease with 
increasing For ^ < 1, the ground state is always characterized by r?, = 0, / = so that 
no vortex is nucleated. For Q = uj, the energy spectrum is made of Landau levels that are 
degenerate in angular momentum. 

The current density in a given eigenstate is defined by 

h , „.m 



{K,i^^n,i - *n,^Vv&;, - 2^-An\^^^') (12) 



2mz^ h 
Its radial component vanishes while its azimuthal component is 



m r h 

We define the angular momentum dependent radius r; by 



> = -(--?Slr)|4'„,,r (13) 



For a given angular momentum state, je{r) is positive for r < ri and hence it gives a 
paramagnetic contribution. It is negative and diamagnetic for r > ri and it vanishes at 
r = r/. 
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B. Spectrum with Dirichlet boundary conditions (DBC) 

We consider now the problem of a trapped rotating boson in a disc of radius R imposing 
Dirichlet boundary conditions. The correspondin g sp ectrum is derived in the same way as 
for the electron in a perpendicular magnetic field jS; 

The Dirichlet boundary condition (DBC) \l/„,i(r = R) = 0, is according to (fTUI) iFi(a, |/| + 
1; $) = where $ = b^^R^ is equivalent to the magnetic flux. The energy spectrum is 
obtained from the zeroes of the confluent hypergeometric function iFi, that can be separated 
into two classes according to the sign of the corresponding angular momentum. For / > 0, 
the energy levels are given by 

^' -^^^ + i;^)=o (15) 

with e = ^"^'^^''^ = e$. For / < the energy levels are given by 

ii^i(^^^^-^,-/ + i;$) = o (16) 

It is important to notice that the equation iFi{a,c]x) = has solutions only for a < 
and, in the interval —p < a < —p + 1, it has exactly p real solutions. This eliminates the 
possibility of having for rotating bosons, a ground state given by the lowest Landau level 
solution for which ^ = 2^^^ ' 

The infinite plane solutions are reached asymptotically for large ^ (see FigO])- There 
is no difference in this spectrum between bulk and edge states. For ^ < 1, the state 
(n, /) = (0,0) is always the ground state (FigCJ). The energy levels become degenerate in 
angular momentum when ^ = 1 and $ ^ 1 (Fig. |21and 3^). Therefore at zero temperature, 
no vortex can be nucleated at - < 1 under the DBC. 

UJ 

We have plotted the energy as a function of the angular momentum for a fixed $ in Fig. 
El For the Landau problem of a charged particle in a transverse magnetic field a similar 
plot, but in a different geometry and only for positive angular momentum states, has been 
used in j^] in order to describe the role of the edge states in the quantum Hall transport. In 
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FIG. 1: Energy spectrum of rotating bosons with Dirichlet boundary condition. The continuous 
lines ( thin lines for n = and thick lines for n = 1) correspond to the energy spectrum derived 
from ij,5)) . Each curve corresponds to a given value of the angular momentum. They correspond 
to n = and / = 0, 1, 2, 3, 4, 5, 6, 7 as well as n = 1 and I = —1,0, 1, 2, 3. The dotted lines represent 
the corresponding infinite plane solutions given by Here is taken to be .75uj. 

the present problem as well as in the Landau problem there is no sharp difference between 
edge and bulk states under DBC. 

IV. CHIRAL BOUNDARY CONDITIONS 

Dirichlet boundary conditions do not provide a way to separate edge from bulk exci- 
tations. As already pointed out, such a separation is necessary for the description of the 
nucleation of a vortex. We therefore propose a set of non-local and chiral boundary con- 
ditions which are more suitable for the present problem. These boundary conditions are 
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FIG. 2: Comparison between the energy spectrum of an electron in a transverse magnetic field 
(dotted lines) and a rotating boson under DBC (continuous lines). The latter is obtained by setting 
^ = 1 in For the curves corresponding to the rotating boson £ is taken to be 0.75. Only 

the lowest Landau level (n = 0) is displayed for both cases. We have plotted the energy levels for I 
between and 7. 



akin to the one introduced by Atiyah, Patodi and Singer (APS) in their study of Index 
theorems for Dirac operators with boundaries j^l- Similar boundary conditions have also 
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42| | . These boundary 



been applied to the Landau problem on manifolds with boundaries 
conditions split the Hilbert space into a direct sum of two orthogonal, infinite dimensional 
spaces with positive and negative chirality on the boundary. The chirality is determined by 
the direction of the azimuthal velocity projected on the boundary. A vortex is nucleated as 
a result of the transfer of a state from the edge to the bulk. 

We have already noticed that for a given value of the angular momentum /, the current 
flows with a different chirality in regions separated by a ring of radius r/ = \/^. We 
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FIG. 3: Edge and bulk states with DBC. The energy levels that correspond to the first two Landau 
levels are shown. The lower level corresponds to n = and the upper one corresponds to n = 1. $ 
and the ratio ^ are respectively taken to 20 and 0.75. 

define the bulk and the edge regions using this particular value of r; as a reference for a 
given angular momentum. The current associated to that particular angular momentum is 
respectively paramagnetic and diamagnetic in the bulk and at the edge. Alternatively, we 
define the bulk and at the edge states for a disc of size R so that bulk states have angular 
momentum / < buR^ whereas edge states have / > b^R^. 

The azimuthal velocity frr^, projected on the boundary of the disc has eigenvalues given 



where $n = b^R^. The chiral boundary conditions are now defined in the following way: 



by 




(17) 
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For A > 0, namely for < < I 



dMR = (18) 

(For any n and henceforth we shall drop the subscript n in -0.) 
For A < 0, namely for / < 

For the first set of wavefunctions that accounts for the edge states we use Neumann 
boundary conditions. We could have used as well Dirichlet boundary conditions. However 
unlike Neumann boundary conditions they give an unphysical discontinuity {4^. These 
wavefunctions are more and more localized towards the outer side of the system for an 
increasing rotation frequency. For states / < that corresponds to wavefunctions localized 
well inside the disc, we impose the mixed boundary condition ()19|) . It is this separation in 
the dispersion relations of the edge states under the choice of CBC that mimics the effect 
of interactions. In the many-body theory of vortex nucleation this se par ation is achieved 
by solving the Bogoliubov equations under various approximations Isfll. 0] and then by 
identifying the collective excitations localized at the surface of the condensate. 



A. Spectrum with CBC 

According to the chiral boundary conditions when $ is increased at a fixed ^, the sign 
of the eigenvalues A(-R) changes from the positive to the negative. Correspondingly the 
energy e of a state with a given n and / changes. This change in energy describes how the 
corresponding state is transferred from the edge Hilbert space to the bulk Hilbert space at 
the point = I ()18|) . This is shown in FigEJ For large $, the infinite plane solutions (|TT|) 
are reached asymptotically. 

For an infinite system we have ^ = n + ^(1 + (1 — ^)0- Therefore for a given ^, is 
a linear function of / with slope (1 — ^). When chiral boundary conditions are applied, this 
behaviour is approximately obeyed for the bulk states. But for the edge states the energy 
increases non-linearly with increasing angular momentum. This is shown in FigO This 
is also the case for Dirichlet boundary conditions (Figl^l). There is however a quantitative 
difference and moreover the bulk and the edge states are now separated under these chiral 
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FIG. 4: Energy levels with chiral boundary condition. Energy levels for the first few angular 
momentum states are shown for n = and n = 1. Each curve corresponds to a given value of the 
angular momentum. They correspond to n = and I between and 5 (thin lines for bulk states and 
thin dotted lines for edge states) as well as n = 1 and I between —1 and 3 with I = —1 corresponds 
to the lowest curve (thick lines for bulk states and thick dotted lines for the edge states). For I > 
the spectrum has a kink at the point I = <I>n ~ 

boundary conditions. For a given n the spectrum is continuous at the point / = but its 
derivative is not. This has been shown in Fig|31for the first two values of n. The derivation 
of the spectrum under these boundary conditions is provided in detail in the appendix. Here 
we focus on the nucleation of vortices under these boundary conditions. 

The Figini shows the effect of an increase of the rotational frequency on the spectrum 
under the choice of CBC. We have plotted the energies of the bulk and the edge states 
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i2/co=0.75 




10 20 

FIG. 5: Edge and bulk states with chiral boundary conditions. Energy levels that correspond to 
n = (lower curve) and n = 1 (upper curve) are shown. Here the edge states are denoted by 
continuous lines whereas the bulk states are denoted by dotted lines. The slope of the energy levels 
is discontinuous at I = where is taken to be 15. 

for four different values of the ratio ^. For each value, the quantity is increased by 
unit steps from 1 to 3 and the corresponding bulk and edge energies are shown. Under 
these conditions, the slope of the bulk energy levels increases while the slope of the edge 
energy levels goes down. The opposite behaviour is observed when, for a fixed the 
ratio ^ increases. Therefore with increasing rotational frequency, states with higher angular 
momenta are transferred from the edge Hilbert space to the bulk Hilbert space. The l-th 
angular momentum state is nucleated in the bulk from the edge by changing from / to 
/ + !. 
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n/co=0.25 
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£1/01=0.5 
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(d) Landau Problem 
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FIG. 6: Effect of a faster rotation. In these figures we have plotted the energies as a function 
of I for a set of ^ values (given above each figure). The three plots in each figure correspond to 
= 1, 2, 3 (the thinnest one for = 1 and the thickest for <I>q = 3). The dotted part corresponds 
to edge states while the continuous part corresponds to bulk states. For ^ = 1 bulk states for all 
three values of fall on the same line. 



The plots that appear in Fig El represent solutions of the stationary Schrodinger equation 
under the choice of CBC at different values of Q and ^^i. To understand the nucleation of 
vortices in the condensate with these solutions we use the fact that a condensate is rotated 
only by nucleating a vortex. Therefore in between nucleations of successive vortices, the 
radius of the condensate remains constant. This is in contrast to the rotation of a rigid 
body which flattens out continuously while increasing the rotational frequency. Let us 
denote by Qi, the characteristic frequency for nucleation of the first vortex. For a given 

< fii, the radius of the condensate is R and it corresponds to $0 = 1. The corresponding 
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FIG. 7: Determination of the characteristic frequency of vortex nucleation. In these figures we have 
plotted the energy ^ against the angular momentum I. The two curves in each figure correspond 
to <I>f7 = 1 (thin), = 2(thick). Corresponding ratios ^ are mentioned above each figure. In each 
curve the continuous part corresponds to bulk states and the dashed-line part to edge states. At 
^ = 0.35 the edge spectrum for <I>q = 1 intersects the corresponding edge spectrum for = 2 (bl 
and b2). At ^ = 0.36, the intersection takes place between the bulk spectrum for = 2 and the 
edge spectrum for = 1 (al and a2). The regions near points of intersection are enlarged in the 
figures a2 and 62. 

bulk and edge energy levels for $n = 1 are obtained from (jl8p and (jl9|) . The condensate 
is therefore the bulk region that contains only the / = state. At = fii, the first vortex 
is nucleated by transferring the / = 1 state from the edge to the bulk so that $0 changes 
from 1 to 2. Therefore the condensate is now defined as the bulk region of $n = 2 and it 
has a larger radius. The rotational flux transferred by nucleating a vortex is thus ^ since 
this is the unit of At f2 = the edge states for $n = 1 intersects the bulk states for 
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$n = 2. This is energetically favoured since a.t Q = Qi, the energy ^(/) of any state with 
/ > 2 is less if $n = 2 rather than = 1. In the corresponding many body description 

0, 

vortices are nucleated when at higher rotational frequency, a surface energy 
barrier disappears. The reduction of the edge state energy with the nucleation of a vortex 
is qualitatively similar to the disappearance of the surface energy barrier. 

The characteristic rotational frequency for the nucleation of the first vortex determined 
using chiral boundary conditions ()18I19|) is between Q = 0.35 u and 0.36 u as shown in 
Fig|7| where cu is the transverse trap frequency. We note that this characteristic nucleation 
frequency is close to the value Q = 0.29 u that has been observed in one of the experiments 
on vortex nucleation All these features cannot be derived from either the infinite plane 
or from DBG. 

To establish a parallel with the case of superconductors let us mention similar results 
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46,43. However 



obtained for the vortex nucleation in a mesoscopic superconducting disc 
in a neutral superfluid there is no Maxwell- Ampere equation that can relate the current to 
the induced rotation . 



V. VORTEX NUCLEATION FOR > J^i 

Thus far we have discussed the nucleation of the first vortex. Here we discuss the bulk 
and edge spectrum under CBC for f2 > We consider the following two cases: 

1. For a further increase of Q beyond Qi, the edge states corresponding to = I crosses 
the bulk states of = / + 1 for / > 2 (FigEl). The corresponding rotational frequencies 
Qi are the nucleation frequencies of the l-th vortices. Each of these vortices is nucleated 
by changing $n to <I>f2 + 1 and it carries one unit of angular momentum. Characteristic 
rotational frequencies are presented in the Table lU Since the bulk and edge energy 
spectra under CBC is obtained by solving the stationary Schrodinger equation in a co- 
rotating frame, solutions at different values of Q can be related one to the other if the 
.oU.o„ . =w.c.ed on ^J.aba^caH. Sue. expen.ent= Q W been perfo..ed 
in a deformed rotating trap where the cross-section of the rotating condensate is an 
ellipse in the plane of rotation. Here we solve the Schrodinger equation in a circular 
domain. This makes the comparison with the experimental results rather difficult. The 

used an approach different from ours, to explain the nucleation 



theoretical work 
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Vortex 


n 


1st 


0.35 - 0.36 


2nd 


0.46-0.47 


3rd 


0.52-0.53 


4th 


0.569-0.57 


5th 


0.603 - 0.604 


6th 


0.629 - 0.630 



TABLE I: Nucleation frequencies of successive (1 = 1) vortices. 

mechanism of the first vortex. It should also be noticed that the adiabatic nucleation of 
successive vortices has not yet been observed experimentally. For a recent theoretical 



work on this issue see |48[ . One of the limitations of our theoretical frame- work is that 
we are unable to determine the spatial arrangement of more than one vortex inside 
the condensate. 

2. Another interesting point is to find the rotational frequencies at which edge states of 
$f7 = 1 intersects the bulk state / = 2 for $n = 3, 4, ■ ■ ■ . The bulk and edge energy lev- 
els of the stationary Schrodinger equations at these rotational frequencies correspond 
to a situation where the rotation is suddenly switched on at a frequency higher than Qi 
and subsequently the system is brought to equilibrium in the co-rotating frame. More 
than one angular momentum states (and rotational fluxes) are transferred to the bulk 
at these characteristic frequencies. This leads to the nucleation of a vortex of angular 
momentum Z > 1 or multiple / = 1 vortices at a time. Within the present theoretical 
framework it is not possible to identify which of these two alternatives is energetically 
favourable. In Table |n] we have listed these rotational frequencies and the number of 
angular momenta transferred to the bulk Hilbert space. Around — = 0.5, the num- 

nn 

ber of vortices nucleated in this way increases very fast. Experimentally [10, uJ| it 
has been observed that under sudden switch-on of the rotation at higher rotational 
frequencies more than one vortices are nucleated and the number of nucleated vor- 
tices is the largest, (see for example Fig. 3 in [ll|) at f2 = This is associated 
with a resonant quadrupolar excitation of angular momentum 1 = 2. The energy of 
this quadrupolar surface mode in units of the transverse trap frequency is determined 
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No. of 


n 


states transferred 




to the bulk 




2 (1,2) 


0.44 - 0.45 


3 (1,2,3) 


0.48-0.49 


4 (1,2,3,4) 


0.494-0.495 


5 (1,2,3,4,5) 


0.498-0.499 



TABLE II: Nucleation frequencies of multiple vortices. In parentheses we mention the angular 
momentum states which are transferred from the edge to the hulk at these values ofQ. 

within the classical hydrodynamic approximation that is very different from our ap- 
proach. Therefore although our findings qualitatively agree with these experimental 
observations, a quantitative comparison with both the experiments and the related 
theory is rather involved. A reason for that is the difference in the trap geometry in 
the two-dimensional plane which is elliptical in the experiment while circular in our 
case. Another reason is the different theoretical approach that have been used to study 
.rfa. e....o. ,„ .Ke .^.-bod. ..eo.y and ou. TKe d,ffe.„ce 

between the characteristic rotational frequencies in the above two cases (Table IJ and 
ITT]) thus shows that the nucleation frequency of a vortex in the condensate depends on 
the number of vortices already present. 



The radius of the condensate is changed when a vortex is nucleated. Since a boson in the 
condensate has one more unit of angular momentum when a vortex is nucleated, its orbit 
has a larger radius. This explains why size of the condensate region increases with successive 
nucleations. When the l-th vortex with one unit of angular momentum is nucleated, the 
condensate region corresponds to $n = / + 1. Therefore at each the radius of the 
condensate region is defined as 



rc\Q=Q, 



'H{1 + 1) 
mVti 



(20) 



Between the nucleations of the /-th and (/ + l)-th vortex, the radius of the bulk region for 

^^^'^^^ and the corresponding edge region grows in size. In Fig|Hl 



decreases as r 



Q u ^^^^^^^^^ , \/ 

using the set of frequencies listed in Table H] we have plotted the condensate radius. In 
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FIG. 8: Size of the condensate as a function of the rotational frequency. The continuous lines gives 
the radius of the condensate r^ as a function of ^ . The dots give the radius of the hulk region for 
<&n = I o,nd for Vii <n < Q-i-i- The length is measured in units of where b^j = 



the same figure we have simuhaneously plotted the radius of the bulk region for $n = /, 
in between the nucleation of the l-th and (/ — l)-th vortices. This plot involves a set of 
disconnected lines and the corresponding differences with the condensate region show how 
the edge region grows between the nucleation of two successive vortices. With the increase 
of the number of vortices in the condensate these two plots get closer one to the other. In 
;he limit of a large number of vortices, the average rotation approaches the rigid body value 



lOl |22| . For example, it has been shown in [10| that experimentally this happens when 
the number of vortices is around 10 jlo| . 

Khawaja et. al. j^l have studied the surface excitations of a trapped three-dimensional 
Bose-Einstein condensed gas using the Gross-Pitaevskii equation and they have associated 
the kinetic energy of the surface region to an effective surface tension. Although our approach 
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is different and that we consider a strictly two dimensional system, we can also define an 
energy associated to the edge by summing the energies of the edge states for $n = /. For a 
finite system then a part of this edge energy is used to nucleate vortices . 

VI. CONCLUSION 

We have presented a one-particle effective theory that describes the process of vortex 
nucleation in a two dimensional rotating condensate using a set of non local and chiral 
boundary conditions. We have shown that the edge states have a dispersion relation distinct 
from the bulk states. These boundary conditions can be understood as a way to mimic the 
effect of the interaction between the bosons. We have demonstrated that the properties 
of vortices thus nucleated agree qualitatively and quantitatively with experimental findings 
although a more thorough comparison is not possible in the absence of a precise mapping 
between this set of boundary conditions and the effective interaction. The expression of 
the characteristic rotational frequency at which vortices are nucleated under adiabatic and 
sudden switch-on of the rotation is calculated and the change in the condensate size with 
increasing rotation is also emphasized. 
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Appendix I 



Case 1: For A>0^0<$n< / plugging the explicit form of the wavefunction (jiup 
in equation (fTHj) we get js^ 

- l)iFi(a, / + 1; $) + ^ iFi(a + 1, / + 2; $) = (21) 
Case 2: For A < there are two possibilities 



Case 2a: For / < ( negative angular momentum) the spectral boundary condition 
(fT^ implies 

+ 4^ + ^^^)^nAr=R = (22) 

or rod 

Using the explicit form of wavefunction given in ()10|1 we find that this demands 

21 iFiia, $) - $(1 - -)iFi(a, |/| + 1; $) = (23) 

UJ 

Case 2b: The other situation is where < / < $n. Again using the explicit form of 
the wavefunction (/ > 0) it can be found that 

iFi(a + 1, / + 2; $) - (1 - -)iFi(a, / + 1; $) = (24) 

It can be seen that the lowest Landau level solution exists only ii Q = uj. If we set 
2bi^ = b,uj = Q where b = and Uc is the cyclotron frequency, all the above results 
matches with the results obtained for the Landau problem on a finite disk under chiral 



boundary condition 



4l|. We solve all these equation with Mathematica j4 



To obtain the characteristic frequencies for the nucleation of successive vortices (Table H} 
we have solved the equation (PT|) for = I and = I + 1 simultaneously. For I = 1 the 
solution yields 0.35uj <Q < .36uj (Fig. Ej). The characteristic frequencies listed in Table ITTl 
are obtained from simultaneously solving equation ()21|) with = 1 and equation with 
<ff7 = 3,4,5,---. 
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